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Abstract 

The system of hydrodynamic-type equations is derived from Alexeev's general- 
ized Boltzmann kinetic equation by two-side distribution function for a stratified 
gas in gravity field. It is applied to a problem of ultrasound propagation and at- 
tenuation. The linearized version of the obtained system is studied and compared 
with the Navier-Stokes one at arbitrary Knudsen numbers. The problem of a gen- 
eration by a moving plane in a rarefied gas is explored and used as a test while 
compared with experiment. It is good agreement between predicted propagation 
speed, attenuation factor and experimental results for a wide range of Knudsen 
numbers 



Introduction 

Fluid mechanics equations in its' most popular form (Euler, Navier-Stokes, Fourier- 
Kirchhoff, Burnett, etc) appear in methods based on Boltzmann kinetic equation by means 
of expansion in Knudsen number (Kn). The first important version of such theory was 
made by Hilbert. It implies analiticity in Kn of both distribution function as well as 
momenta functions. Further development of the theory by Chapman-Enskog and Grad 
jl] weaken analiticity condition of the momenta on Kn. It allowed to deride NS, Burnett 
equations for the Chapman-Enskog method and 13 momenta Grad equations widely used 
in fluid dynamics description. Failures in deep Knudsen regime penetration recovered 
by direct attempts with many-moment theories lead to more deep understanding of the 
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problem |2l El E] • The Knudsen independent expansion of the basic (Boltzmann) equa- 
tion was used, namely one of Gross- Jackson, starting from the celebrating BGK model. 
The unification of Chapman- Enskog and Gross- Jackson approaches j5] exploits an idea of 
nonsingular perturbation method in its Frechet expansion form jH]. 

One of important verification of fluid dynamics system relates to the problem of sound 
propagation. Its simplest version considers the plane harmonic wave with the correspon- 
dent dispersion relation. Such case obtained by linearization of the basic system repro- 
duces the known experiments of [7] rather well. It incorporated in a direct scheme of 
kinetic approach [Hj. The fluid mechanics systems, based on BGK jH] model of collision 
integral, obtained recently in |H] and jTHj, give good results for velocity of sound in Kn 
0.1 -7- 10 but fail in attenuation description [TP. 

Developing the method based on Gross- Jackson collision integral for a non-isotropic 
fluid for a problem, which specifies J2] a direction in it we use an idea of von Karman 
to divide the phase speed with respect of particle velocity direction along/against the 
direction axis ^^l- Such situation takes place if a gas is stratified in gravity field, that 
yields appearance of interne gravity waves branch with the obvious necessity to account 
wide range of Kn J3] . 

Struchtrup ^Sl regularizes 13-moment Grad equations doing the same thing as a 
test. His linearization results in a dispersion relation, which acoustic branch gives an 
attenuation coefficient that also does not fit experiments. 

Our article is devoted to this problem; we tried to improve the results on a way of 
next Gross- Jackson model [16j, the tendency was good but the changes were not enough. 
Considering an alternative possibility to compensate the discrepancy in relaxation timees- 
timation, we adress to the Alexeev generalization of Boltzmann equation ^7] . 

Alexeev-Boltzmann equation looks like: 



locity and radius vector of the particle, respectively, r is the mean time between collisions, 
is the collision Boltzmann integral. 

We apply our method for the generalized Boltzmann equation of Alexeev and such 
"joint" theory gives a better agreement with the experimental data [7] for attenuation at 
arbitrary Knudsen number. 

Generalized fluid dynamics equations 

Consider the kinetic equation with the model integral of collisions in BGK form 0: 
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- local- equilibrium distribution function, = y2kT/m denotes the average thermal 
velocity of particles of gas, u = z/q exp(— - is the effective frequency of collisions 
between particles of the gas at height z, H = kT/mg - is a parameter of the gas stratifi- 
cation. It is supposed, that density of the gas n, its average speed U = {Ux, Uy, Uz) and 
temperature T are functions of time and coordinates. 

Following the idea of the method of piecewise continuous distribution functions let's 
search for the solution / of the equations(^ as a combination of two locally equilibrium 
distribution functions, each of which gives the contribution in its own area of velocities 
space: 



(3) 

f- = n- — — exp ' , Vz<0 



here n^, U^, depending on t, z are functional parameters. 

The double number of parameters of the distribution function results in its deviations 
from a local-equilibrium one. In the range of small Knudsen numbers Kn << 1 we 
should have = n~ ,T~^ = T~ , = U~ and distribution function start from a local 
equilibrium and at the small difference between the functional 'up' and 'down' parameters 
produces the Navier-Stokes equations. The theory is also valid at big Kn(free molecular 
regime) [T8]. 

We restrict ourselves by the case of one- dimensional disturbances U = (0, 0, U), using 
a set of linearly independent momenta functions: 



(4) 



Here the first three functions are collisional invariants. Let's define a scalar product 
in velocity space: 

< (fn, f > = < (fn >= J dv (pnf ■ (5) 

<m>= p , <mVz>= pU , < imC^ >= |^/cT , 

where ^ = V — U is the peculiar velocity. Here p = nm is mass density, Pzz is the diagonal 
component of the pressure tensor, Qz is a vertical component of a heat flux vector, Qz is a 
parameter having dimension of the heat flux. 

If we now multiply the kinetic equation with the model integral of collisions in BGK 
form by ipi and integrate over velocity space, the fluid dynamic equations appear: 
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-r— fp[/^ + lOU'^P,, + 20U^q, + 5 < m^,^ > f/+ < m^,^ > 

= -2uq,-3uU{P,,~-^kT), 

m 

where 

^3 = f < e.' >, J4 = f < efe^ > 

The system (jT)) of the equations according to the derivation scheme is vahd at all Kn. 
To close the description it is enough to plug the two-side distribution function into 
that yields for n^, U^, as function of p, U, T, Pzz, qz, Qz- We base here on an expansion 

in small Mach numbers M = maxl — I , up to the first order. The values of integrals (jHl) 

Vt 
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Figure 1: The inverse non-dimensional phase velocity as a function of the inverse Knudsen 
number. The results of this paper are compared to Navier-Stokes, Chen-Spiegel p!U] . 
regularization of Grad's method ^Sl and the experimental data [7] 



as functions of thermodynamic parameters of the system (|7j) are: 

■h = --p{ ) + -. + -{—) pToT 
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m m rrv- „ rrr' „ 

kTo kTo m m 

J4 = Q—qz + ^—qz + QpoU^ - 6pU^. 

Substitute Jj into ((Tj) gives modification of fluid dynamics equations at arbitrary 
Knudsen numbers. 

Let us report the results of investigation of quasi-plane waves parameters as a function 
of Kn. For this purpose we proceed in a standart way: linearizing the system ((7j) that 
impose the dispersion relation as a link between frequency and complex wave number. 

In the model of hard spheres in the continual limit r can be connected with the 
dynamical viscosity rj 2J 

TP = 0.786?7. 

r and u are linked : 

In figures 1, 2 a comparison of our results of numerical calculation of dimensionless 
sound speed and attenuation factor depending on r ~ is carried out in a parallel way 
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Figure 2: The attenuation factor of the hnear disturbance as a function of the inverse 
Knudsen number. 

with the results by other authors. The Navier-Stokes prediction is quahtatively wrong 
at big Knudsen number. Our results for phase speed give the good consistency with 
the experiments at all Knudsen numbers. However, our results for the attenuation of 
ultrasound are good (as we can see in experiment) for the number r up to order unity 
and in the free molecular regime. Taking into account disadvantages of model integral of 
collisions it is planed to consider kinetic equation with full integral of collisions. It will 
permit to describe processes in transition regime. 
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